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Understanding the behaviour of topologically ordered lattice systems at finite temperature is a way of assess- 
ing their potential as fault-tolerant quantum memories. We compute the natural extension of the topological 
entanglement entropy for T > 0, namely the subleading correction J topo to the area law for mutual information. 
Its dependence on T can be written, for Abelian Kitaev models, in terms of information-theoretic functions and 
readily identifiable scaling behaviour, from which the interplay between volume, temperature, and topological 
order, can be read. These arguments are extended to non- Abelian quantum double models, and numerical results 
are given for the D(Sz) model, showing qualitative agreement with the Abelian case. 

PACS numbers: 



The notion of topological order was first introduced in the 
context of the fractional quantum Hall effect []]]. It aims at 
identifying phases that cannot be separated by local order pa- 
rameters. Such phases can exhibit exotic phenomena such as 
topologically protected ground state degeneracy or quasipar- 
ticle excitations, called anyons, with statistical behavior that 
is different from bosons or fermions yy. Besides the fact 
that they reveal unusual states of matter, topologically ordered 
systems are also interesting because they might allow for in- 
trinsically fault-tolerant quantum computation |[3L 0] - In such 
systems, the division of a quantum algorithm as initialisation, 
unitary evolution and read-out |5] would translate into creat- 
ing anyons, braiding them, and fusing them back to the vac- 
uum respectively. While conceptually appealing, the robust- 
ness of this topological quantum computation against realistic 
noise models is far from being fully assessed. 

This paper is devoted to quantifying how temperature af- 
fects a topologically ordered medium. For that, we will use 
an entropic topological order parameter, / t opo, and focus on 
lattice spin systems that are exactly solvable ^ lfl6ll . We will 
show that, at any fixed temperature, J topo is nonzero only if 
the size of the system is finite. Importantly, we will exhibit 
a scaling relation describing how a given increase of the sys- 
tem size can be compensated by a vanishing decrease of the 
temperature. After recalling some notions on the topological 
entropy, we will provide a complete analysis for the simple 
toric code model. Then, we will turn to more general mod- 
els, show how to compute entropic quantities, and provide nu- 
merical evidence for scaling laws in the case of the simplest 
non-abelian quantum double model, the one based on D(Ss) 
Jit]. We believe that our findings are relevant to topological 
quantum computation |01 ■ 

We start our discussion by briefly recalling the notion of 
topological entropy |7, 8]. It will play a central role in the 
following. Let us consider a two-dimensional many-body 
quantum system in a pure state of its ground subspace and 
let R : R c denote a bipartition of this system, with R be- 
ing connected. Let us further assume that our system satis- 
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Figure 1: Division of a torus (or a sphere) into four regions. 



fies an 'area' law. That is, the von Neumann entropy of re- 
gion R grows linearly with the size of its boundary, Sr = 
—trpR In pr = a'\dR\ — 7' + S(R), where pr denotes the re- 
duced density matrix corresponding to region R, a 1 is a con- 
stant, and where lim^R^^ S(R) = 0. As first pointed out 
in Jal, and further analysed in 0,01, a non-zero value for the 
constant correction 7' reveals topological order. Indeed, 7' 
is a topological invariant of the system and is related to the 
quantum dimension T> of the model at hand as 7' = In V. 

Let us consider a system defined on a torus or on a sphere, 
and divide it into four regions A, B, C, D (see Fig[T|). We will 
use 

hoyo = 1a + Ib + Ic — Iab — Iac — Ibc + Iabc- (1) 

as a topological order parameter. (Ir = Sr + Sr? — Sr u r^ 
denotes the mutual information between a region R and its 
complement R c .) This definition of 7 topo amounts to replace 
von Neumann entropies by mutual informations in the defini- 
tion of topological entropy introduced in |8y. At zero temper- 
ature, itopo = —27'. Our choice is motivated by the fact that, 
at finite temperature, the mutual information between a region 
R is a measure of the correlations of this region with its envi- 
ronment. The lattice systems we are going to study obey an 
area law: 

I R = a\dR\- 1 , (2) 

and have a finite correlation length £ jjj. Therefore, 

a quantity where correlations due to a finite £ cancel out, and 
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that reveals correlations due to topological order only, as the 
topological entropy, 7', at zero temperature 10]. Direct substi- 
tution of Eq. (12]) into Eq.fHJ shows that 7 topo = 7 lfl7ll . 

We now want to compute I t opo at finite temperature for the 
toric code yfl. Consider a torus T tiled into L x L square pla- 
quettes. Let us associate a two-level system ('spin') with each 
edge of this tiling and let us assume that these spins interact 
through the hamiltonian 



H 



(3) 



where the index p (resp. s) runs over all plaquettes (resp. ver- 
tices) of the tiling. The operators B p and A s are defined as 
Bp = Yl iep <Ji and A s — Y\ ies crf- The hamiltonian H is 
made of local terms, all commuting with each other and, as it 
turns out, its eigenstates can have arbitrary values for the op- 
erators B p and A s , up to the constraint J\ B p = Y[ s A a = 1. 
Therefore, these eigenstates are labelled by a triple of quan- 
tum numbers: \<j>, c, w). A pattern <f> denotes the position of 
all plaquette or "flux-type" excitations, a pattern c indicates 
the position of all vertex or "charge-type" excitations, and w 
indexes the topological degeneracy of the state for a fixed con- 
figuration of defects. The latter quantum number is made of 
two bits (wi and W2) that label the values of the integrals of 
motion of 2 -operators around non-contractible loops on the 
torus (Wilson loops) [3] . We have 



H\<f>,c,w) = (E + 2\<f>\ + 2|c|) \<f>,cw), 



(4) 



where Eq = — 2L 2 is the ground state energy, and \<p\ 
(resp. |c|) denotes the number of flux excitations (resp. 
charge excitations) of the pattern <b (resp. c). If Pi de- 
notes the projector onto the sector of energy Ei = Eq + 4i 
and di = trPj denotes its dimension, we have that di = 

^X/n^,,n c <L 2 /2 Sn^+n e =i (2714,) (2n c ) ' 

We now compute all the thermodynamic al quantities we 
need, assuming that our system is at thermal equilibrium at 
inverse temperature (3. Some details of our calculations will 
be presented elsewhere IToll . The partition function reads 

Z(fi,L) = tre"^ = Eflo e ~ 0Eid i = ((2 cosh/3) L2 + 
(2sinh/3) L ) 2 , and the von Neumann entropy of the state 
p = e~@ H /Z(/3, L) of the whole torus is derived from the 
identity 



S m = lnZ-(P/Z)dZ/d/3. 



(5) 



In order to compute the von Neumann entropy of a region 
R, and eventually I topo , we first observe that pr(4>,c,w) = 
tTflJ'/'jC, w){4>,c, w\ (R c = T\R) depends on w only if R 
contains non-contractible loops (w cannot be measured lo- 
cally). Also, up to total charge conservation, pn(<f>,c,w) 
will only depend on the excitations located on plaquettes and 
crosses inside R. Therefore, we distinguish between three 
types of flux excitations: those located on plaquettes with sup- 
port on R, cf)R, those located on plaquettes with support on 
Rc, <Pr c , an d those located on plaquettes with support par- 
tially on R and partially on R c , 4>dR- Similarly, we divide the 



charge excitations into cr, cr c , cqr. The crucial observation 
that helps to compute pr = tiR c e~P H /Z(f3, L) and eventu- 
ally S(pr) is that <9i?-fluxes and 9i?-charges can be driven 
inside R c by application of Pauli operators acting on links of 
R c , that is: |</>, c, w) is of the form Ur c \4>r, 4>r ,cr, cr c ,w), 
for some unitary operator Ur c . To lighten the notations, let 
us denote by and c\r c the patterns of excitations (both flux 
and charge) in R and in R c respectively. All excitations in 
R (resp. R c ) can be fused into a single excitation q R (resp. 
q R J through a unitary operator U R (qn, q R ) ® U' Rc (q Rc , q R J 
that relates |q#, q^ e , w) and a state \qR,q\,w), such that 
R and R c contain at most one excited plaquette (vertex). 
As in iH, we will call site a plaquette and an adjacent ver- 
tex. The possibly excited plaquette and the possibly ex- 
cited vertex within R (resp. R c ) can always be chosen to 
form a site. The state \q R , q R , w) can be obtained, from a 
ground state |£,u>), by a product of unitary operators along 
a string connecting the site where q R is located to the site 
where q R is located. In summary, |qw) can be written as 
Ur(c[ R ) ® UR c (qR c ,q aR )\£„w) for some unitary operators 
Ur(cir) and U Ra (q flc , q aR ). (q 9i? refer to the charges which 
are neither fully in R, nor fully in R c .) 

We are now in a position to compute Sr. The thermal state 
of the toric code reads 



e -0(E o +AE\q\) 



(6) 



where AE = 2 is the energy associated with a single excita- 
tion, and where q = q^ U qi? c U qoR- Therefore the reduced 
density matrix of the system R reads p R = J2 W Eq H , qSc : q aR 

ftR [UR(ciR)\£,,w)(£,,w\UR{ciRy}. To sim- 



-/3(E +AE|q]) 
Z{P,L) 



plify further this expression, we observe that two reduced 
states pr(4>, c, w) and pr(4>', c', w') are orthogonal whenever 
cr) 7^ c' R ). Indeed two such states can be perfectly 
distinguished through a measurement of A s or B p operators 
having support on R. Therefore, the sum Ylq R t s actually a di- 
rect sum. If we assume that R is contractible in neither of both 
directions on the torus, then the values of the Wilson loops w 
can be revealed by measurements having support fully on R 
and the sum J2 W a ^ so turns to be a direct sum. So, pr, can be 
written as 

Pr = C(q R )tiR c UR(ciR,x)\£,w)(£,w\U R (<iR,x), 

(7) 

„-/3(B +2|q SR |+2|q Rc |) 



where 4C(q/?) 



4 e -2f!<iBl) 



Z(J3,L) ^qas,qH c 

is the marginal probability of a configuration of defects qR. 
It is the decomposition (0 that allows to compute Sr exactly. 
From it, we find that Sr separates into a ground state area 
contribution and a finite j3 contribution 



Sr — St 



V(p > N p (R),L) + V([3 ) N*(R),L) ) (8) 



where S R 



S(tr Re \t,w)(t,w\) = (\9R\ - l)ln2 M, 
and where N p (R) (N*(R)) denotes the number of plaquettes 
(crosses) fully contained in R. The function V can be com- 
puted exactly using elementary combinatorial identities fioll . 
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When pr is fully contractible, the direct sum over w appear- 
ing in (f7]i should be replaced by a simple sum. As a result, the 
expression for Sr picks a — In 4 additive correction. 

We have used Eq.© to compute I t opo- We have found that 
at any finite (3, I topo vanishes in the limit where L — > oo, a 
result that is consistent with those of Ref. lfTTll . and indicates 
that (i) the toric code exhibits no temperature-driven phase 
transition for itopoi an d that (ii) the toric code is likely not a 
scalable quantum memory in the strictest sense JU as far as 
temperature-induced errors are considered lfl8ll . Our compu- 
tations also reveal that the mutual information between a re- 
gion R and its complement satisfies an area law of the form 
Eq.©. Computing 7 topo from it in the limit L — > oo, and 
considering |<9i?| = AvL, with v < 1, we have found a re- 
markably simple expression: 



/topo = 2 In 2 - h 2 (p e ; en ) - h 2 (pl 



(9) 



where ti2(x) = — xln(x) — (1— x) ln(l-x) denotes the Shan- 
non entropy of a binary outcome random variables, p p ma ~ 
p even ~ (]_ + Q" 2 l 2 )(\ + Q^-" 2 ) L ' 2 )/2 ( 1 + L *) is the prob- 
ability that the region R contains an even number of excited 
plaquettes (crosses), with 6 = tanh/3. Eq.© allows to un- 
derstand simply why Ttopo vanishes at finite temperature when 
L — > oo: in this limit, p p ven ~ p e J en ~ 1/2. This equation also 
reveals a scaling law: at fixed value of v, I topo only depends on 

p and L through the parameter t = 6 .In particular, a fixed 
value of t (and thus a fixed value of the topological mutual 
information), corresponds to the following relations between 
size and temperature: 



/3(t,L) =]nL 



&T(t,L) 
dL 



iln(I In i) + 0(ZT 2 ), (10) 



+ 0{L- 3 ). 



L{\aL - i ln(| In j) + 0(L~ 2 )) 2 



(ID 

This last relation tells us how an increase of the size of the 
system should be compensated by a decrease of temperature 
in order to maintain a fixed value of the topological entropy. 
We understand it as an important nuance over the fact that 
/topo vanishes when L — > oo, at finite temperature. In partic- 
ular, it shows that the rate at which the temperature should be 
decreased, in order to maintain a fixed value of 7, decreases 
with the size of the system. 

We now turn to a family of non-abelian models, those based 
on a quantum double. From now on, we will consider systems 
defined on an oriented lattice with the geometry of a sphere. 
A quantum degree of freedom with basis states labelled by 
the elements of a finite group G is associated with each link 
of this lattice. These links interact through a hamiltonian of 
the form (0), with vertex opertators A s and plaquette opera- 
tors B p that still commute. One could write down their form 
explicitly 01, but it will not be necessary here. A natural way 
to deal with such models would be to start by providing a de- 
scription of the complete set of eigenstates of H similar to the 
one we have used for the toric code. However, when G is non- 
abelian, diagonalising H seems to be a difficult problem. Yet, 



we can argue that we actually do not need to as far as we are 
only interested in the topological mutual information. The el- 
ementary excitations of H live on sites (a site is a combination 
of a vertex and an adjacent plaquette). We will restrict to that 
part of the spectrum of H such that excitations are elemen- 
tary and pinned at fixed non-adjacent sites. This restriction 
can be thought of as additional error correction, where some 
plaquettes and vertices are over-protected so that they never 
get excited by thermal fluctuations (or only with vanishing 
probability). Therefore, we believe that 7 t opo for this modified 
model can only be larger than for the full spectrum. 

The space of n excitations pinned at fixed sites has the 
structure ^ 



'In • 



(12) 



where each index ^ runs over all possible quasiparticle types 
for site i. Eq.dTZli simply means that different excitation pat- 
terns lead to orthogonal states. Each space < K qi ... qn splits as 



®0Cq n ®M qi ... qn , 



where the spaces % qi correspond to the local degrees of free- 
dom of the quasiparticles |3J. The fusion space, M. gi ... q „, 
is what makes non-abelian anyonic systems so special. For 
abelian models, for which the fusion rules are trivial, the di- 
mension of the fusion space is equal to one. A non-trivial fu- 
sion space appears when fusing two anyons can yield different 
quasiparticles J5t]: 



<7, 



„ Xq b = Y. N a 



b Qc- 



(13) 



The system has to fulfill some neutrality conditions, i.e. its 
state should be such that fusing all the particles yields the 
trivial particle, denoted 1, with certainty |[l2tl . The dimen- 
sion of M. qi ...q n depends on the tensor N*^ as follows 0]: 



dim M c 



N bl TV, 2 
qiq2 b iq3 



K 



We 



l qi...q n Sb 1 ...6 n _ 2 ' 

observe that computing this quantity amounts to contracting 
a (quasi) translationally invariant matrix product state lfl4ll . 
Likewise dimJf [n] can be computed efficiently, as well as the 
partition function of the model 



Z[fi, n) = ]T . . . J2 II <%)dim M qi ... q 



-PE(qi...q n ) 



'In 



where E{q\ . . . q n ) is the energy associated with a defect con- 
figuration qi . . . q n , and where d(qi) is the dimension of the 
space % qi . 

To compute I to po, we consider a situation where a pair of 
anyons \qq) is created in such a way that anyon q lies in some 
region R and its anti-quasiparticle q lies in the complementary 
region. The von Neumann entropy of region R then reads (H] 



S(p R ) = S(p s ^) + \nd q . 



(14) 



The entropy of a region when the system is in a thermal 
state can be computed once we are able to calculate the en- 
tropy of a region when the system is an arbitrary defect con- 
figuration. In turn, just as for the toric code, the latter entropy 
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reduces to computing the entropy when there are only two 
anyons in the system, and one lies inside the region we are 
interested in. Let n R and n Rc denote the number of sites con- 
tained in region R and R c respectively (n R + n Ro = n). Let 
1r = <?i • • • Qjir label the configurations of types of anyons 
living on the sites contained in region R. The total fusion 
space splits as: 



I 

topo 



-2 




M **,4*. = M L ® K' Rc ® K» (15) 



where M„ R denotes the fusion space associated with event 
where all quasiparticles fuse to an anyon of type b. 

Due to total anyonic charge conservation, dinuM^ r = 1. 
The decomposition (JT3J induces the following representation 
of the thermal state of H: 



Ah = 







y{(l R -> b,H\,s R ;q[ R -> b,fi 2 ,s R J, 



(16) 



where the projector T(q R -» 6, fix, s R ;q! Rc -> b,fi 2 ,s Ra ) 
refers to a pure state where all (resp.q^. ) quasiparticles 
contained in i? (resp. i? c ) fuse to b (resp. 6) through the 
channel fix, 1 < fix < dimM^ R (resp. fi2, 1 < fii < 
dimMt ). £(q fl ) = E?=i denotes the energy asso- 

ciated with the configuration q R , and s R is a collective index 
for the internal degrees of freedom of the quasiparticles con- 
tained in region R, 1 < s R < d(q R ) = Y[j=i d{q 3 )- E{c^ Ra ) 
and s Rc are defined likewise for the region R c . Note that the 
ground state energy is now made equal to zero, upon shifting 
the hamiltonian by a multiple of the identity . 

The von Neumann entropy of the reduced state of the region 
R now reads 



Z R {p,b)Z Ra {P~b) 



In dh 



-kxZ Ro (0,b))]+J2 



Y R (f3,b)Z Rc (P,b) 
Z{(3,n) 



-\xxZ{f3,n), (17) 



where Z R (f3,b) = ^d^e-^^dimM^, and 
where Y R ((3, b) = (3-^Z R ((3, b). The von Neumann entropy 
of the whole s phe re is given by Eq.©, while the ground state 
entropy reads lH : Sf = In \G\{\dR\ - 1). 

Eq.(fT7b has allowed us to study numerically how the topo- 
logical mutual information behaves as a function of j3, for 
G = S3, the smallest non-abelian finite group. Our results 
are shown on Figf2] The systems we have considered are four 
tiled spheres, all with 96 x 96 plaquettes. The first sphere con- 
tains n = 64 sites, the second 144, the third 256 and the fourth 



Figure 2: Topological quantum mutual information as a function of 
the inverse temperature /3 for the D(Sa) model. (Color online.) 



576. Although these systems are small, they are large enough 
to show that non-abelian systems are affected by temperature 
in the same way as the toric code: (i) For a fixed number of 
sites, there exists a finite width topological phase, (ii) / topo 
tends to for small values of (3, and to —3.58352 = —2 In 6 
for sufficiently large values of (3, as expected since the order 
of S3 is six. (iii) The larger the number of sites, the larger 
the value of j3 where the topological mutual information van- 
ishes. Finally, we have observed that when I lopo is plotted as a 
function of ne~ 2 P, the curves collapse when n becomes large. 
We take this observation as strong evidence that the behaviour 
of /topo for the -D(Sa) model is governed by the same scaling 
variable as the toric code. Indeed, ln(6> L ) ~ L 2 e~ 2f3 in the 
limit where (3 and L are large (n oc L 2 for the toric code). 
We therefore believe that the discussion held after Eq.© also 
holds in this case, and more generally for any quantum double 
model. 

In conclusion, we have shown that the interplay of ther- 
mal effects, lattice size, and topological order (as measured 
in the subleading correction to the area law) is encoded in 
well-defined scaling relations (TTOb - dTTT > . In particular, the rate 
with which the temperature should be decreased to fight the 
effect of thermal fluctuations vanishes in the limit of large 
lattices. These relations seem to hold for both Abelian and 
non-Abelian systems. As a byproduct, we have derived a for- 
mula for the entropy of a region for non-abelian quantum dou- 
ble models defined on a lattice. This formula depends on the 
model only through the fusion tensor N*^ and through the en- 
ergy associated with each quasiparticle. It is therefore tempt- 
ing to believe that it holds in a more general context. Finally, 
it is appealing to try to give an operational meaning to 7 to po, 
by connecting its value with the possibility of using a topolog- 
ically ordered medium as a robust quantum memory. But to 
the best of our knowledge, such a connection is still an open 
problem, even at zero temperature. Eq.© might contribute 
establishing it. 
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